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1 ■ ABSTRACT 

; I ' Using extended sets of cosmological hydrodynamical simulations of galaxy clusters, 

we present a detailed study of scaling relations between total cluster mass and three 
^ ' mass proxies based on X-ray observable quantities: temperature of the intra-cluster 

medium, gas mass and the product of the two, Yx = MgasT. Our analysis is based on 
two sets of high-resolution hydrodynamical simulations performed with the TreePM- 
SPH GADGET code. The first set includes about 140 cluster with masses above 5 x 
10"'^^ h~^MQ, with 30 of such clusters having mass above 10^^ h~^AlQ. All such clusters 
I have been simulated in two flavours, with both non-radiative physics and including 

' cooling, star formation, chemical enrichment and the effect of supernova feedback 

^SJ , triggering galactic ejecta. The large statistics offered by this set of simulated clusters 

is used to quantify the robustness of the scaling relations between mass proxies and 
total mass, to determine their redshift evolution, and to calibrate their intrinsic scatter 
and its distribution. Furthermore, we use a smaller set of clusters, including 18 halos 
with masses above 5 x 10^'^ h~^MQ, four of which are more massive than 10^^ h~^MQ, 
to test the robustness of mass proxies against changing the physical processes that 
are included in the simulations to describe the evolution of the intra-cluster medium. 
Each cluster is simulated in seven different flavours to study the effect (i) thermal 
}_j , conduction, ('ii^ artificial viscosity, (^nij cooling and star formation, (^wj galactic winds, 

03 . and (v) AGN feedback. 

As a general result, we find the M-lx scaling relation to be the least sensitive 
one to variations of the ICM physics, its slope and redshift evolution being always 
very close to the predictions of the self-similar model. As for the scatter around the 
best-fitting relations, its distribution is always close to a log-normal one. Mgas is the 
mass proxy having the smallest scatter in mass, with values of a\nM — 0.04-0.06, 
depending on the physics included in the simulation, and with a mild dependence on 
redshift. As for the mass-temperature relation, it is the one with the largest scatter, 
with (7\nM<^ 0.1 at z = 0, increasing to ^ 0.15 at z = 1. The intrinsic scatter in the 
M-Yx relation is slightly larger than in the M-Mgas relation, with a\nM — 0.06 at 
z = and 0.08 at z = 1. These results confirm that both Yx and Mgas mass proxies 
are well suited for cosmological applications of future large X-ray surveys. As a word 
of caution, we point out that the analysis presented in this paper does not include 
observational effects expected when measuring temperature from fitting X-ray spectra 
and gas mass from X-ray surface brightness profiles. A detailed assessment of such 
effects will be the subject of a forthcoming paper. 

Key words: Cosmology: Miscellaneous - Methods: Numerical - X-rays: Galaxies: 

^ . Clusters. 
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1 INTRODUCTION 

Measuring the mass of galaxy clusters represents the most 
important aspect for the use of such structures as tracers 
of cosmic evolution. The cosmological parameters depend, 
indeed, on the evolution of cluster's mass function a nd on 
large-sc ale clustering properties (e .g., 'Borga ni et al.l 



ter mass re.g..lEvrard et al.lll99d iBartelmann fc Steinmeta 



"Schuecker et al.l |2002| : IVoid l2005al : Vikhlini n et al 



2001 



200C 



Mantz et al.ll201Cl[ ). Restricting the discussion to X-ray ob- 



servations, to estimate cluster masses we need to assume 
that the intra-cluster medium (ICM) sits in hydrostatic 
equilibrium within the cluster potential well, generally as- 
sumed to have spherical simmetry. Further, we need to ob- 
tain gas density and temperature profiles from high-quality 
X-ray data. Limits on this approach are represented not only 
by the strong assumptions made but also by the difficulty 
of obtaining data of sufficient quality. This is especially true 
for distant clusters, at zf^ 0.5, for which observations of the 
needed sensitivity would be prohibitively expensive in terms 
of required exposure time. 

To overcome this problem one generally resorts to the 
observational determinations of mass proxies, i.e. of suitably 
observable quantities easier to measure and related to total 
cluster mass through the so-called "scaling relations". Ide- 
ally, a robust mass proxy should be characterized by a low 
intrinsic scatter with cluster mass, implying that it neeeds to 
be relatively insensitive to the cluster dynamical state. Dif- 
ferent examples of X-ray mass proxies have been used over 
the last decade to derive cosmological constraints from X- 
ray clusters surveys. The easiest-to-measure mass proxy is 
the total X-ray luminosity, Lx. While this quantity is known 
to co r relate with cluste r mass (e.g.. iReiprich k, Bohringeij 
l2002l : iPratt et all bOOgl ). the Lx-M relation has a fairly 
large scatter, of about 40 per cent, due to the sensitivity 
of X-ray luminosity on the details of the gas distribution 
in central regions and, therefore, on the dynamical state of 
the cluster. Other commonly adopted cho ices of mass prox- 
ies a r e the X-ray temperature Tx ( e.g. |Henrv fc Arnaud 



199ll : iMarkevitchI 1 19981 : llkebe et al. I l2ooa r iPierpaoli et al 



20031 ) and the cluster gas mass M |;aB(e.g. IVikhlinin et al 



20031 : IVoevodkin fc VikhlininI |2004| ). In order to estimate 
these quantities detailed X-ray brightess profiles must be ob- 
tained, requiring longer exposure times than those required 
by Lx only. However, they are expected to be more robustly 
linked to total collapsed mass. 

In this context cosmological hydrodynamical simula- 
tions are playing an increasingly relevant role in calibrat- 
ing the mass proxies, in understanding the systematics on 
cluster mass measurements, and in defining new mass prox- 
ies. Simulations provide the most advanced theoretical tool 
to capture the complexity of the hierarchical build-up of 
galaxy clusters and of th e gas-dynamical processes (e.g. 
iBorgani fc Kravtsovl [ioogl . for a recent review on cluster 
simulations). A typical example is provided by the role 
that simulations played in the last years to calibrate the 
expected level of violation of hydrostatic equilibrium in- 



iRasia et al.l |2004|: 


Nagai et al.l l2007l: iJeltema et al. 


lAmeglio et al.l 200£ 


: Piffaretti et al.ll2010l: Lau et al. 



20091 ). Another example is given by the application of 
simulations in the study of the accuracy of X-ray tem- 
perature and of gas mass as tracers of the total clus- 



Il99e ; iMuanwong et aklbood : lEttori et al]|2004l : iNagai et al 
i2007r ). With the increasing numerical resolution and degree 
of realism of the physical processes included, these calibra- 
tions are becoming trustworthy as convincingly shown by 
different works. A good match between simulations and ob- 
serva tions was shown for temperature profiles of clusters 
(e.g..lLoken et al.ll2002l:lBorgani et al.ll2004l:[Kav et al.ll2007l : 
iNagai et all bOO?!: iLeccardi fc Molendil |2008|). g as de nsity 



profiles (e.g., Roncarelli et al.l 2006 Croston et al.llioo^ ') and 



pressure den sitv fe.g.. I Arnaud et ak l bood). 

These comparison are mostly restricted to the region 
within -Rsoo Q, after excluding the central regions affected 
by the complex physical processes regulating the "cool core" 
structure of the ICM. 

As for the cluster mass proxies simulations predict that 
the Tx-M relation has an intrinsic scatter which is quite sen- 
sitive to the presence of subs tructures in the ICM and to the 
clust er dynamical state (e.g.. lO'Hara et al]|2006l : lYang et al.l 
l2009ll . The gas mass proxy is instead expected to be less sen- 
sitive to the cluster dynamical state, thus implying a smaller 
scatter in its scaling relation with total mass. Moreover it 
has the advantage that it can be essentially measured from 
X-ray imaging alone. There is however a caveat, since the 
dependence of gas mass fraction on cluster mass, cluster- 
centric distance and on redshift are still not yet completely 
und erstood (e.g., Ettori et"all bood ) . 

iKravtsov et al.l (|20oi r introduced the X-ray equivalent 
of the integrated SZ flux, the Yx parameter, defined as the 
product of the gas mass with the cluster temperature. From 
a set of sixteen simulated galaxy clusters these authors found 
it to be a low-scatter (5-8 per cent) mass proxy, due to the 
anticorrelation of the residuals in the scaling relation of tem- 
perature and gas mass with total mass. Furthermore, the 
scaling of Fx with total mass was found to be in close agree- 
ment with the prediction of the simple self-similar model, 
based on the assumption that gas follows dark matter and 
gravitational effects only determine the thermal content of 
the ICM (e.g.. K aiser ,1986 ). The Yx-M relation is rela- 
tively insensitive to clu ster mergers , a re sult that was in- 
dependently verified bv lPoole et al.l (|2007l ) and by Rasia et 
al. (2010) with simulations of cluster mergers. A number 
of recent observational works estimated Yx and its correla- 
tion with X -ray luminosity. Data from bot h XMM-Newton 
( Pratt et a l. 2009) and Chandra telescope (jVikhlinin et al.l 
I2OO9I : lMaughan.2007. ) showed that the observed Yx-M re- 
lation has slope and redshift evolution in agreement with 
predictions of the self-similar model. With their sample of 
both nearby and distant clusters out to z ~ 1, they confirm 
that gravitational processes are indeed responsible for the 
bulk of the total thermal content of galaxy clusters. 

A more recent analysis of scaling relatio n s for simulated 
galaxy clusters presented by IStanek et al.l (|2010l ) did not 

confirm the anti-correlation of the residuals of M Mgas 

and M—Tmw They found that the lowest-scatte r mass proxy 
is Mg as- A similar result was also found by lOkabe et al.l 
(|2010l) who computed the scaling relation between X-ray 



^ We define Ma as the mass contained within the radius R/^ en- 
compassing an overdensity of A times the cosmic critical density, 
Pciz) = 3H{zf/8TTG. 
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mass proxies and the total cluster mass obtained from lens- 
ing observations of 12 clusters at a ~ 0.2. From their analysis 
the conclusion is that Mgas has the lowest intrinsic scatter 
of ~ 10%, with respect to the ~ 20% of the Yk proxy. 

At present, Tx, Mgas and Ix are considered the most 
robust indicators of cluster mass and, in fact, the most recent 
cosmological applications of X-ray cluster surve ys have been 
based on the use of such mass proxies (e.g. , Ilfenrv et al.l 
l2009l : IVikhlinin et al.|[2009l : iMantz et al.ll2010l '). 

Ideally, rather than relying on simulations, it would be 
preferrable to calibrate mass proxies directly from obser- 
vational data. An observational calibration, while possible 
in principle, requires highly precise measurements of cluster 
mass in a way which is completely independent of X-ray ob- 
servations, e.g. from gravitational lensing or from the study 
of the cluster internal dynamics as traced by member galax- 
ies. Indeed, any determination of the scatter between any of 
the above three X-ray mass proxies and cluster mass based 
on the application of hydrostatic equilibrium would be just 
inconclusive, owing to the fact that the relation is between 
two highly correlated quantities. On the other hand, precise 
measurements of lensing masses are no w available only for a 
relatively small number of objects (e.g . , iMahdavi et al.|[2008l : 
IZhang et allboiol : [Okabe et al.ll201Gl ). Although progresses 
in this direction are expected in the future, simulations rep- 
resent today a valid alternative to calibrate mass proxies 
and to understand their robustness. 

A careful study of mass proxies through simulations re- 
quires a large enough sample of clusters simulated both with 
adequate resolution and by including a range of physical 
processes. The most direct way to achieve large statistics of 
clusters would be to carry out simulations of large cosmolog- 
ical boxes (e^g. 'Borgani et al."2004'; 'Gottloeber et al."2006f; 
iKav et al.1 1200 7: Burns ct al. 2008; Hartley et al. 2008). 
However, these simulations have to severely compromise be- 
tween box size, achievable resolution and details with which 
different physical processes (e.g., star formation and feed- 
back from different sources) can be studied. As an alter- 
native, one could simulate at high resolution only specific 
regions surrounding clusters, which are previously identified 
from low-resolution simulations of large cosmological vol- 
umes. The advantage of this approch is that it can provide 
a more realistic numerical description, in terms of resolu- 
tion reached and of physical processes included, only within 
the "zoomed-in" Lagrangian regions surrounding clusters. 
As a matter of fact, this resimulation procedure is non- 
trivial to implement, thus limiting the statistics of resimu- 
lated cluster s so-fa. r presented in the literature. For instance, 
iDolag et al.l l|2009l ) presented a set of 18 such clusters, of 
which however onl y 4 have M 2 00 > 10^^ h ~^M(^ (see also 
iFabian et aLll2010l ). iLau et al.1 (j2009l ) and IPuchwein et al.1 
l|2008h presented sets of resimulated clusters of comparable 
size, again with only a couple of them being as massive as 

In this paper we use two sets of simulated clusters, 
whose combination provides us with both a large statis- 
tics of massive objects and a range of physical processes 
included over which to test the stability of the scaling re- 
lations. The first set of clusters (Set 1 hereafter) contains 
about 140 clusters, 30 of which have masses larger than 
10^^/i~^Mq. For all of them, we carried out simulations 
with both non-radiative physics and by including cooling. 



star formation and galactic winds powered by Type-II super- 
nova (Sn-II) explosions. The second set of simulations (Set 2 
hereafter) includes 18 clusters, with 4 of them hav i ng ma sses 
larger than lO^'^ h-^M@. We refer to lDolag etlH (|2009l ) for 
a complete description of the second set of simulations. This 
smaller set of simulated clusters complements the first set 
since it has been simulated by changing in seven different 
ways the description of physical processes (e.g. artificial vis- 
cosity, cooling, star formation, feedback efficiency from SN 
and from gas accretion onto supermassive black holes) . The 
complementary analysis of these two cluster sets allows us to 
address in detail the stability and robustness in the calibra- 
tion of X-ray mass proxies with simulations, by combining 
statistics of massive halos, numerical resolution and range 
of physical processes included. 

In this first paper we will carry out an analysis of scaling 
relation of mass proxies against mass by neglecting all ob- 
servational effects. For this reason ICM temperature will be 
mass-weighted (where not listed differently) , and all quanti- 
ties will be estimated in a three-dimensional analysis with- 
out including projection effects. In this way, we will assess 
the intrinsic performances of Tx, Mgas and Ix as tracers of 
the true mass, while we defer to a forthcoming paper (Rasia 
et al., in prep.) a detailed analysis of the impact of observa- 
tional biases. 

This paper is organized as follows. We present in Sec- 
tion [2] the simulations analysed in this pap er. After a de - 
scription of the GADGET simulation code (|Springell [20051 ) 
and of the physical processes included in the simulations, 
we describe the initial conditions and the general charac- 
teristics of the resulting samples of simulated clusters. In 
Section [3] we present the results of our analysis. After in- 
troducing the definitions of mass proxies, we present results 
at 2 = 0, along with an assessment of the dependence of 
the scaling relations on the physical processes included in 
the simulations. We finally discuss the evolution of scaling 
relations and of their intrinsic scatter with redshift. We dis- 
cuss our results and summarize the main conclusions of our 
analysis in Section |4l 



2 SIMULATIONS 

2.1 The simulation code and adopted physical 
models 

In this Section we introduce the simulation code and all the 
different physical processes considered in re-simulating our 
objects. Simulations have been ca rried out using the Tree- 
PM SPH GADGET l|Springell2005l ') code. Oldest GADGET 2 
ve rsion was us e d for Set 2 simulations, already presented 
bv lDolag etai] (|2009l '). while Set 1 was completed using the 
newest version 3. The difference between the two versions re- 
sides mainly in the different algorithm adopted for domain 
decomposition. Both versions of the code use segments of 
the space-filling Peano-Hilbert curve to decide the particles 
to be assigned to different processors. Unlike GADGET 2, 
the newest GADGET 3 version allows each processor to be 
assigned also disjoined segments of the Peano-Hilbert curve. 
This turns into a substantial improvement of the work-load 
balance assigned to the different processors, especially for 
simulations, like those presented here, in which the compu- 
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tation cost is largely concentrated within a quite small frac- 
tion of the physical volume of the computational domain. 

We add here below a brief description for each of the 
physical models adopted within our simulations. 

NR-SV (non-radiative, standard viscosity): 
non-radiative runs, using the standard reference scheme for 
artificial viscosity implemented in GADGET This prescrip- 
tion is based on th e formulation of a rtificial viscosity orig- 
inally presented bv iMonagh an' ("1997*), also inclu ding a vis- 
cosity limiter as proposed by Balsara (Il995l ) and ISteinmet j 
(|l996h . 

NR-RV (non radiative, reduced viscosity): 
non-radiative runs, using r educed time-dependent visc osity, 
as originally proposed by iMorri s fc Monaghan' and 
implemented in GADGET by [Oplag ct al. (2005). In this 
scheme, artificial viscosity decays away from shock regions, 
thus allowing a higher degree of turbulence in the velocity 
field to develop. 

CSF (cooling star formation and feedback): 
runs including radiative cooling for a zero-metallicity 
plasma, including heating/coo ling from a spatially un iform 
and evolving UV background (|Haardt fc Madaull 19961 ). Gas 
particles above a given threshold density are treated as mul- 
tiphase, so as to provide a sub-resolution description of the 
inter-stellar medium, accord ing to the model described by 
ISpringel fc Hernauisd (|2003h . Within each multiphase gas 
particle, a cold and a hot-phase cohexist in pressure equi- 
librium, with the cold phase providing the reservoir of star 
formation. Kinetic feedback is implemented by giving some 
SNII energy to gas particles in the form of kinetic energy, 
thus mimicking galactic ejecta powered by SN explosions. In 
these runs, galactic winds have a mass upload proportional 
to the local star-formation rate. We use ~ 340kms~^ 
for the wind velocity, which corresponds to assuming half of 
energy released by SN-II being co nverted to kinetic energy 
for a Salpeter IMF l|Salpeteij|l955l ). 

CSF-C (csf and thermal conduction): 
the same as the CSF runs, but also including the e f fect o f 
ther mal conducti o n, as described by lJubelgas et al.l (|2004l ) 
and iDolag et al.l (|2004l ). The conduction efficiency is as- 
sumed to have a value of 1/3 of the Spitzer conductivity. 

CSF-M-W (csf, metals and galactic winds): 
radiative runs with star formation described through the 
same multiphase model of the CSF runs. In addition, CSF-M- 
W runs also include a description of metal production from 
chemical enrichment contributed by S N-II, SN-Ia and AGB 
stars, as described bv lTornatore et aO (2007). Stars of differ- 
ent m ass, distributed according to a Salpotor IMF (Saltaeterl 
[l95l), release metals over the time-scale determined by 
the corre sponding mass-dependent life-times (taken from 
IPadovani & Mattcucci 1993). The metallicity dependence of 
radiative cooling is included by using the cooling tables by 
[Sutherland fc Dopital l|l993D . In these runs, the velocity of 
galactic ejecta is assumed to be v^, = 500 km s"'^, thus cor- 
responding to a more efficient SN feedback that in the CSF 
runs. 

CSF-M-NW (csf, metals and no winds): 
the same as CSF-M-W but with no winds, i.e. by excluding 
the effect of kinetic feedback from galactic winds. 

CSF-M-AGN (csf, metals and AGN): 
the same as the CSF-M-W runs, but replacing SN feedback in 



the form of galactic winds, with the effect of AGN feedback 
triggered by gas accretion onto supermassive black holes 
(BH). Detail s abou t this feedback model are discussed in 
iFabian et all l|2010l ). while we summarize here only the main 
points. The sc heme is a slig h t mo dification of the model 
introduced by ISpringel et al. l|2005l ) , where BH are repre- 
sented by sink particles initially seeded in resolved DM ha- 
los. These particles increase their mass by gas accretion and 
merging with other BH particles. Eddington-limited Bondi 
accretion produces a radiated energy which correspond to a 
fraction = 0.1 of the rest-mass energy of the accreted gas. 
A fraction e/ of this radiated energy is thermally coupled 
to the surrounding gas. We use £/ = 0.05 for this feedback 
efficiency, which increases to e/ = 0.2 whenever accretion 
enters in the quiescent "radio" mode and takes place at a 
rate smaller than one- hundred of the Eddi ngton limit (e.g. 
ISiiacki fc Springellbood : iFabian et al.ll201Gl ). 

We remark here that galaxy clusters in Set 2 were sim- 
ulated with all the presented baryon physics models, while 
simulations of Set 1 clusters were performed with NR-SV 
and CSF-M-W physics only. 

The diversity and heterogeneity of the description of 
physical processes determining the evolution of the ICM is 
not meant to provide a systematic study of the parameter 
space describing the ICM physics. However, quantifying the 
variation of the scaling relations between mass proxies and 
true cluster mass for each of such models will allow us to as- 
sess the robustness of mass proxies against the uncertainties 
in the description of the ICM physics. 

2.2 Initial conditions and simulation sets 

Our analysis is based on two sets of simulations of galaxy 
clusters. The two sets, while having similar mass and force 
resolution, are quite different in a number of aspects and 
allow us to perform complementary tests on the robustness 
of mass proxies. The large statistics of Set 1 clusters allows 
us to calibrate with precision the slope and scatter of the 
scaling relations along with their evolution. While the vari- 
ety of different physical processes of Set 2 are used to assess 
the robustness of mass proxies against the uncertainties in 
the description of the ICM physics. 

We provide in the following a more technical description 
of these two sets of simulations. 

2.2.1 Set 1 

is based on simulations of 29 Lagrangian regions extracted 
around as many clusters identified within a low-resolution 
N-body cosmological simulations. We provide here below a 
short description of these initial conditions, while a more de- 
tailed presentation will be provided in a forthcoming paper 
(Bonafede et al., in preparation). The parent DM simula- 
tion follows 1024'^ DM particles within of a box having a co- 
moving side of lh~^Gpc. The cosmological model assumed 
is a fiat ACDM one, with Q.m = 0.24 for the matter den- 
sity parameter, fibar = 0.04 for the contribution of baryons. 
Ho — 72kms~^Mpc~^ for the present-day Hubble constant. 
Us — 0.96 for the primordial spectral index and ag = 0.8 
for the normalization of the power spectrum in terms of the 
r.m.s. fluctuation level at z = within a top-hat sphere of 



© 0000 RAS, MNRAS 000, 000-000 



Mass proxies from simulations 5 



radius. With the above choice of parameters, this 
cosmological mo del is consi stent with the CMB WMAP-7 
constraints ([Komatsu et alj 2010). The selected Lagrangian 
regions have been chosen so that 24 of them are centered 
around the 24 most massive clusters found in the cosmo- 
logical volume, all having viriaU mass 

Further five regions have been chosen around clusters in the 
mass range Mvir ~ (1 - 7) x 10^* fe'^M©. 

Within each Lagrangian region we increased mass res- 
olution and added the relevant high-frequency modes of 
the power spectrum, following the Zoom ed Initial Condi- 
tion (ZIC) technique l|Tormen et al.|[T997^ . Initial conditions 
for hydrodynamic simulations have been generated by split- 
ting each particle within the high-resolution region into a 
DM and a gas particle, having mass ratio such to reproduce 
the cosmic baryon fraction. Outside high-resolution regions 
resolution is progressively degraded, so as to save compu- 
tational time, while preserving a correct description of the 
large-scale tidal field. We generated initial conditions at dif- 
ferent, progressively increasing, resolution. For the hydrody- 
namic simulations presented here the mass of each DM parti- 
cle in the high-resolution region is moivi — 8.47x10* h-^Me, 
with mgaa — 1.53 x 10* h~^MQ for the initial mass of gas 
particles. Using an iterative procedure, we have shaped each 
high-resolution Lagrangian region in such a way that no 
low-resolution particle contaminates the central "zoomed 
in" halo at z = at leeist out to 5 virial radii of the cen- 
tral cluster. This implies that each region is sufficiently large 
to contain more than one interesting cluster with no "con- 
taminants" out to at least one virial radius. In total, we 
find ~ 140 clusters with M500 ^ 5 x lO^/i'^M©, out of 
which about 40 have Afvir > 5 x 10^" H'^^Mq and 30 have 
Mvir > lO^'' h'^^MQ. We show in Fig.[T]the cumulative mass 
distribution for this set of simulated clusters. 

Simulations have been carried out using a Plummer- 
equivalent softening length for the computation of the grav- 
itational force in the high-resolution region which is fixed to 
e = 5 /i~^kpc in physical units at redshift z <2, while being 
kept fixed in comoving units at higher redshift. As for the 
computation of hydrodynamic forces, we assume the SPH 
smoothing length to reach a minimum allowed value of 0.5e. 

Clusters from this set have been simulated using the 
NR-SV and CSF-M-W models. The NR-SV set is the only 
one simulated with GADGET 3 code. 



2.2.2 Set 2 

is based on 9 Lagrangian regions selected from a par- 
ent DM only c osmo logical box with a size of 479/i"^Mpc 
lYoshida et al] l|200ll ). The cosmological model assumed is 
ACDM with SIm = 0.3 for the matter density, h = 0.7 for 
the Hubble parameter, f2bar = 0.039 for the baryon density 
parameter and ff8=0 .9 for the normaliz ation of the power 
spectrum. We refer to lDolag et al.l (|2009l ) for a more detailed 
description of this set of clusters. Briefiy, this set includes 
overall 18 clusters with Afvir > 5 x 10^^ h"^ Mq, out of which 



Here we define the virial mass Mvir, as the mass contained 
within the radius encompassing the overdensity of virialization, a s 
predicted by the spherical collapse model (e.g.. lEke et al.lll996l) . 
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Figure 1. The cumulative distributions of virial masses for clus- 
ters in Set 1 (red solid histogram) and in Set 2 (blue dashed 
histogram). We only count "clean" clusters, for which no con- 
taminant low-resolution DM particles are found within their cor- 
responding virial radii. 



only 4 have Mvir > lO^'^ H'^Mq. The mass distribution of 
these clusters is shown in Fig. [T] 

Mass resolution is quite close to that of Set 1 (mDM — 
1.9 x 10* /i~^M0 and rrigas ^ 2.8 x 10^ /i~^A/q for DM and 
gas particles respectively), with also the same choices for the 
size of the gravitational softening scale and for the minimum 
allowed value of the SPH smoothing length. Clusters from 
this set have been simulated with GADGET 2 using all the 
seven above described different models for the description of 
the baryon physics. 



3 RESULTS AND DISCUSSION 
3.1 Mass proxies 

The mass proxies that we consider in our analysis are the 
mass- weighted temperature Tmw , the gas mass Mgas and the 
product of gas mass and temperature, Yx ~ Mgas x Tmw. 
All these quantities are computed within the fiducial ra- 
dius -R500 , which typically corresponds to the outermost ra- 
dius where scaling rel ations from observational data analy- 
sis are provided (e.g., Ivik hlinin et a"l] |2009l : ISun et al.|[2009l : 
iPratt et al.ll2009l : iMaugh an 2007). With the purpose of re- 
producing the procedure of observational analyses we ex- 
clude the central regions, R < 0.15J?5oo, in the computation 
of the temperature. The reason for this choice in the analy- 
sis of real clusters is that temperature profiles show a rather 
high degree of diversity in such centr al regions depending 
on th e ir degree of "cool-coreness" (e.g., ^Leccardi fc Molendil 
I2OO8I : IVikhlinin et~aLl I2OO6I : IPratt et al., ,2007. ). Excluding 
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Figure 2. Scaling relations at 2; = between the total mass computed within -R500 and three mass proxies: Mgas, Tmw and Yx. The 
evolution of the relation is scaled accordingly, with 7 = 0, 1 and 2/5, respectively. On each figure we plot the results for non-radiative 
(blue open squares) and radiative (red filled circles) runs for clusters in Set 1 with M^^^ > 5 X 10^^ h~^MQ. The short dashed (red) 
and dashed (blue) lines are the best fit relations of respectively radiative and non-radiative cluster simulations with slopes fixed to the 
self-similar value. Left panel: scaling with Mgas; central panel: scaling with mass— weighted temperature, Tmw; right panel: scaling with 
Yx • The values of Tmw appearing in the central panel and entering in the computation of Yx in the right panel are computed within the 
radial range (0.15— l)-R5oo • 
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Figure 3. Distributions of the residuals in log— mass in the scaling relations at 2 = between mass proxies and Mtot,500i for the non- 
radiative (upper red histograms) and radiative (lower blue histograms) runs of simlated clusters from Set 1. Residuals are computed 
with respect to the relation that fits at best the simulated clusters data. Left panel: scaling with Mga^; central panel: scaling with 
mass-weighted temperature, Tmw; right panel: scaling with Yx- Dotted curves in green show the Gaussian distribution having the same 
variance of the distribution of residuals in mass (as reported in Tables [T] |2] and |3)l . 



such regions suppresses the scatter in the scal inR relation 
involving temperature (e.g., IPratt et al.l |2009| . and refer- 
ences therein). Furthermore, it places us in the regime where 
simulations are in c loser agreement with observations (e.g., 
iBorgani et al.ll2004l ). The choice of the radial interval where 
to compute mass proxies represents our only attempt to in- 
troduce "observational effects" in our analysis. Indeed, in 
this paper (Paper 1) we prefer to analyse in detail the in- 
trinsic performances of the three above mass proxies, i.e. 
neglecting any possible observational bias which will be ad- 
dressed in Paper II (Rasia et al., in prep.). We want to stress 
out that in this analysis i) we estimate all quantities in three- 
dimensions (i.e., neglecting projection effects); ii) unless oth- 
erwise stated, we use the mass-weighted definition of tem- 
perature; iii) we do not remove gas clumps at relatively low 



temperature when we compute both temperature and gas 
mass. 

The simplest description of scaling relations between 
galaxy cluster mass and X-ray observables is that p rovided 
by th e self-similar model originally proposed by l|Kaiseij 
Il986l ). According to this model, the thermodynamical prop- 
erties of the ICM are determined only by the action of grav- 
ity, under the a ssumption o f virial equilibrium and spherical 
simmetry (e.g.. lVol3l2005bl . for a review). Within this model 
cluster total mass is the only parameter that determines all 
the thermodynamical properties of the hot intra-cluster gas. 
Accordingly, the cluster gas mass within the radius Ra is 
proportional to the total mass computed within the same 
radius, 

Mtot,A = CMgMgas,A , (1) 
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with constant of proportionaUty CMg being independent of 
redshift. The total mass is related to the temperature ac- 
cording to 

Mtot,A = E{z)-^CtTI^'' , (2) 

where E{z) = H{z)/Ho defines the evolution of the 
Hubble parame t er. T herefore, as defined in the paper by 
iKravtsov et alJ (|2006l ) the self-similar scaling of Yx with 
cluster mass is given by 

Mtot,A = Cy^ E{z)-^/"Y^% . (3) 

In the above relations involving ICM temperature the mass- 
averaged Tmw should be used, since it is more directly related 
to the total ICM thermal content. As such, this is the quan- 
tity relevant for the predictions of the self-similar model, 
rather than the spectroscopic temperature which is more 
affected by the complexity of the ICM thermal structure. 

In the following we analyse the scaling relations for sim- 
ulated clusters with masses M^ir > 5 x 10^^ h~^MQ, at five 
different redshifts, z = 0, 0.25, 0.5, 0.8 and 1. To fit the scal- 
ing relation between total mass, Mtot,500, and a given ob- 
servable mass proxy X500, we use a power law 

M,„,5oo = Cx(^)" . (4) 

Following the work bv lKravtsov et ahl ()2006l l the normaliza- 
tion factor Xo was fixed to 2 x lO^M©, 3 keV and 4 x lO" 
keV Mq for Mgas, Tmw and Yx, respectively. We also per- 
form the fits by fixing the slope ax to the self-similar value 
of ass = 1-0, 1.5 and 0.6 for M500 - Mgas, M500 — Trnw and 
M500 — Yx, respectively. 

The values of the slope ax and normalization Cx are 
obtained by performing a best-fitting on the log-log relation 
by using the least -squares Marquardt-Levenberg algorithm 
l|Press et al.lll992l ). This algorithm calculates at each itera- 
tion the sum of the squared differences with a new set of pa- 
rameter values. The Marquardt-Levenberg algorithm selects 
the parameter values for the next iteration. The process con- 
tinues until a preset criterium is met, either the fit has con- 
verged (the relative change in the residuals is less than 10~^) 
or it reaches a preset iteration count limit. The intrinsic scat- 
ter in natural logarithm of the total mass, InMtot, around 
each best-fit scaling relation is computed as the quadratic 
difference between the single cluster data and the relation 
fitted on the whole cluster set. 



3.2 Scaling relations at z = 

In this section we will discuss the behavior of M500 — Mgas, 
M500 — Tmw and M500 — relations for the radiative (CSF- 
M-W) and non-radiative simulations (NR-SV) at 2 = for 
the simulated clusters of Set 1. 

In Fig. [2] we show the relations between the total clus- 
ter mass M500 and Mgas (left panel), Tmw (middle panel), 
Yx (right panel). Besides the results of the analysis for ra- 
diative (CSF-M-W in red circles) and non-radiative (NR-SV 
in blue squares) simulations, we plot the best-fit relations 
obtained when the slope is fixed to the corresponding self- 
similar values, shown with the red long-dashed and the blue 
short-dashed lines for the CSF-M-W and NR-SV simulations, 
respectively. 



Galaxy clusters in non-radiative simulations show a be- 
haviour close to the self-similar one for the three studied 
relations. Indeed, when we consider the slope as free param- 
eter the new best fits (not shown) over all NR-SV clusters 
return values consistent with ass, namely, a = 0.981±0.004 
for Mgas, a = 1.517 ± 0.012 for Tmw and a = 0.597 ± 0.03 
Yx, respectively. Results of the fitting relations are repoted 
in Tables [TJ [2I and |3] for the Yx, Mg^s, and Tmw, respec- 
tively. The agree ment with the self-s imilar scaling for Tmw 
was found also bv lStanek et~all (|2010l ) for clusters extracted 
from the non-radiative version of the Millenium Gas Simu- 
lations (MGS). 

When radiative cooling and star formation are included 
in CSF-M-W simulations, the average quantities change thus 
modifying slope and normalization of the scaling relations. 
The CSF-M-W clusters (red circles in Fig. [2]) deviate from 
the self-similar scaling relation for both Mgas and Tmw. As 
for the M500 — Mgas scaling, the best fit relation has a shal- 
lower slope. This is due to the fact that for radiative sim- 
ulations the conversion of baryons into stars is relatively 
more efficient in lower mass clusters, thus corresponding to 
a lo wer mass fraction of gas left in the hot diffuse phase 
fe.g.lBorgani et al.ll20oi : IKravtsov et al.ll2005l : iFabian et all 
|201(J ). As for the M500 — Tmw relation, the results for the 
radiative runs are characterized by a slope steeper than the 
self-similar one, with a — 1.615 ± 0.016. In this case, the 
effect of including radiative physics is that of increasing the 
temperature of the ICM by a larger amount in less massive 
systems. Indeed, the more efficient cooling in lower-mass 
systems causes a relatively stronger adiabatic compression 
of gas in central core regions, as a consequence of the lack 
of pressure support. As a result, temperature profiles in ra- 
diative simulations of clusters arc relatively steeper in less 
massive systems (see also Borgani ct al. 200j). 

The deviations with respect to self-similarity that we 
detect in the M500 — Mgas and M500 — Tmw relations for ra- 
diative simulations are not present in the M500 — Yx rela- 
tion. In fact the best fit slope of cluster data is 0.591 ±0.003, 
close to the value of ass = 0.6 predicted by the self-similar 
model. As apparent from the right panel of Fig. [S] the effect 
of including cooling, star formation and SN feedback in the 
CSF-M-W runs has only a modest effect of the M500 — Yx 
scaling. The reason for this is that the Yx mass proxy pro- 
vides a measure of the total thermal content of the ICM. 
Since this total thermal content is dominated by the mech- 
anism of gravitational accretion of baryons within the DM- 
dominated potential wells of clusters, it is not surprising that 
the scaling of Yx with total cluster mass very close to the 
self-similar prediction and that its value, at a fixed cluster 
mass, is weakly sensitive on the inclusion of non-radiative 
physics. This is even more true when we exclude central 
cluster regions, as we did in our analysis, which are more 
affected by radiative cooling. 

Besides the stability of the scaling relations against vari- 
ations of the ICM physics, another useful quantity to judge 
the robustness of mass proxies is the intrinsic scatter around 
the best fitting relations. In fact, the intrinsic scatter and 
its distribution are also an important quantities to look at 
when deciding which mass proxy is the b est one for cosmo- 
logical applicatio ns of galaxy clusters (e.g. [Lima fc Hull2005l : 
'Sha w et al.ll201 (f). Looking at Fig. [2] we note that the mass- 
weighted temperature is the proxy with the largest scatter. 
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2 = z = 0.25 z = 0.50 z = 0.80 z = l 

NR-SV best— fitting parameters 

ay^ 0.597 ± 0.003 0.599 ± 0.003 0.604 ± 0.003 0.601 ± 0.004 0.604 ± 0.005 

logCy^ 14.190 ± 0.002 14.188 ± 0.002 14.189 ± 0.002 14.177 ± 0.003 14.185 ± 0.004 

o-inM 0.064 0.067 0.067 0.070 0.073 

Tir^Yx 0-107 0.111 0.111 0.117 0.122 

NR-SV self— similar scaling (ay^ = 0.6) 

logCy^ 14.191 ± 0.002 14.188 ± 0.002 14.187 ± 0.002 14.177 ± 0.002 14.183 ± 0.003 

o-i„M 0.064 0.067 0.067 0.070 0.074 

crir^Yx 0.107 0.111 0.112 0.117 0.123 

CSF-M-W best-fitting parameters 

(XYx 0.591 ± 0.003 0.590 ± 0.002 0.596 ± 0.003 0.591 ± 0.003 0.596 ± 0.005 

logCy^ 14.249 ± 0.002 14.247 ± 0.002 14.250 ± 0.002 14.248 ± 0.003 14.252 ± 0.004 

CTinM 0.050 0.051 0.054 0.061 0.076 

o-lnYx 0.084 0.086 0.091 0.104 0.128 

CSF-M-W self-similar scaling {cty^ = 0.6) 

logCvx 14.250 ± 0.003 14.250 ± 0.002 14.252 ± 0.002 14.253 ± 0.002 14.255 ± 0.002 

o-inM 0.052 0.054 0.055 0.063 0.076 

(Tin Fx 0.087 0.089 0.091 0.104 0.127 

Table 1. Fitting parameters for the M500 — relation for simulated clusters of Set 1. Results are obtained through a log-log linear 
regression usmg eq. g}, with Fx.o = 4 X lO^^keVM©. Reported errors on the fitting parameters correspond to Icr standard deviation in 

the linear fit. Also reported are the values of the r.m.s. scatter in Mtot.soo and in Yx,500 around the fitting relations. Different columns 
refer to results at different redshifts. First two lines report results for non-radiative (NR-SV) simulations, the other two instead the 
results for radiative (CSF-M-W) simulations. For each simulation we report first the best-fitting parameters and below the results when 
the slope of the scaling relation is fixed at the value predicted by the self-similar model. 

2 = z = 0.25 2 = 0.50 2 = 0.80 2 = 1 

NR-SV best— fitting parameters 

"Mj,„, 0.981 ± 0.004 0.983 ± 0.003 0.991 ± 0.004 0.991 ± 0.004 0.999 ± 0.005 

logCMj,„, 14.154 ± 0.002 14.142 ± 0.002 14.136 ± 0.002 14.130 ± 0.002 14.127 ± 0.002 

(TinM 0.055 0.048 0.047 0.041 0.040 

o-inMga, 0.056 0.049 0.047 0.041 0.040 

NR-SV self-similar scaling (oy^^ = 1) 

logCMg„, 14.155 ± 0.002 14.144 ± 0.002 14.138 ± 0.002 14.132 ± 0.001 14.128 ± 0.001 

o-inM " 0.061 0.052 0.048 0.042 0.040 

o-inMga, 0.061 0.052 0.048 0.042 0.040 

CSF-M-W best-fitting parameters 

"Mga, 0.929 ± 0.003 0.926 ± 0.003 0.929 ± 0.003 0.924 ± 0.004 0.933 ± 0.005 

log Cmj,„, 14.310 ± 0.002 14.298 ± 0.002 14.291 ± 0.002 14.285 ± 0.002 14.288 ± 0.003 

(TinM 0.039 0.043 0.039 0.046 0.047 

o-inMga, 0.042 0.047 0.042 0.050 0.050 

CSF-M-W self-similar scaling (ay^ = 1) 

log Cm<,„, 14.317 ± 0.004 14.320 ± 0.003 14.320 ± 0.002 14.320 ± 0.002 14.322 ± 0.002 

o-inM 0.090 0.092 0.078 0.080 0.068 

(TinMg^^ 0.090 0.092 0.078 0.080 0.068 

Table 2. The same as in Table [Tl but for the A/soo — Afgas relation. In this case we use Mgas.o = 2 X IO^'^Mq for the normalization 
factor in eq. (|4]l. 
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0.25 



: 0.50 



: 0.80 



NR-SV best— fitting parameters 

"T„„ 1-517 ± 0.012 1.515 ± 0.016 1.534 ± 0.016 1.498 ± 0.021 1.495 ± 0.027 

log Ct„„ 14.513 ± 0.006 14.524 ± 0.007 14.535 ± 0.007 14.507 ± 0.009 14.530 ± 0.010 

crinM 0.120 0.145 0.134 0.146 0.150 

crinT„„ 0.079 0.096 0.087 0.098 0.101 

NR-SV self— similar scaling (ctYx — 1-5) 

log Ct„„ 14.508 ± 0.005 14.519 ± 0.005 14.524 ± 0.004 14.508 ± 0.005 14.531 ± 0.005 

o-inM 0.120 0.145 0.135 0.146 0.150 

o-inT„„ 0.080 0.097 0.090 0.098 0.100 



CSF-M-W best-fitting parameters 

"T„„ 1.615 ± 0.016 1.614 ± 0.013 1.640 ± 0.017 1.608 ± 0.019 1.596 ± 0.028 

log Ct„„ 14.427 ± 0.006 14.436 ± 0.005 14.450 ± 0.006 14.445 ± 0.006 14.444 ± 0.009 

o-inM 0.111 0.106 0.127 0.132 0.163 

o-inT„„ 0.069 0.066 0.077 0.082 0.102 

CSF-M-W self-similar scaling (oy^^ = 1.5) 

log Ct,„„ 14.413 ± 0.007 14.410 ± 0.005 14.414 ± 0.005 14.417 ± 0.004 14.419 ± 0.005 

(7inM '"™ 0.134 0.127 0.147 0.142 0.168 

(TinT„„ 0.089 0.085 0.098 0.094 0.112 



Table 3. The same as in Table [T] but for the Afsoo 
eq.im. 



- Tmw relation. In this case we use Tmw.o = 3 keV for the normalization factor in 



both for CSF-M-W and NR-SV runs. To display graphically 
the scatter of each relation we plot in Fig. [3] the distribu- 
tion of residuals, with respect to the best-fit relation, in 
log-mass. In all panels the distribution of the distribution 
of the residuals is close to a log-normal one. The intrin- 
sic mass scatter in the residuals of the M-Mgas relation is 
ainM ^ 0.039 and 0.055 for the CSF-M-W and NR-SV simu- 
lations, respectively (see also Table [2]). The largest scatter is 
obtained for the mass-weighted temperature Tmm, reaching 
values of crinM — 0.11 and 0.12 for non-radiative and radia- 
tive simulations, re spectively, thus similar to what found by 
IStanek et all (|2010l ) for their non-radiative simulations (see 
also Table [31). As for the Ix proxy, the intrinsic scatter in 
mass is ainM ^ 0.050 and 0.064 for CSF-M-W and NR-SV 
simulations, respectively. 

Therefore, while the Ix mass proxy is found to have 
a scaling with mass which is closest to the self-similar be- 
haviour, almost independently of the inclusion of radiative 
physics and SN feedback, the proxy having the lowest scatter 
ag ainst mass is Mpas . T his result is in line with that obtained 
bv IStanek et al.l ( 2010l ) from the analysis of clusters simu- 
lated with SPH, with both non-radiative physics and includ- 
ing cooling plus the effect of a phenomenological pre-heating. 
Similarly to our analysis, Stanek et al. did not attempt to 
include observational eflects in their analysis, while they 
analysed scaling relations at a larger radius, i?200 instead 
of i?5oo hke in our analysis. We also note that indications 
for a lower scatter in the M500 — Afgas relatio n than in the 
M500 — Fx relation have been also reported bv lOkabe et al.l 
(|2010l ) from the analysis of observational data. Their anal- 
ysis was based on XMM-Newton data for the measurement 
of X-ray mass proxies and on weak lensing data for the mea- 



surement of total cluster masses. Th ese results are however 
at variance with the original result bv lKravtsov et al.l (|2006l ) 
who found instead Fx to be the lowest-scatter mass proxy 
from their analysis of cluster simulations based on the ART 
adaptive-mesh refinement code. In their analysis, Kravtsov 
et al. showed that the residuals AT in the M500 — Tmw rela- 
tion and A Mgas in the M500 — Afgas relation anticorrelate 
for their radiative simulations when computed with respect 
to the best-fitting self-similar relation, thus justifying the 
smaller scatter found in the M-Yx_ scaling. Indeed the latter 
can be computed as ay = o"t„„ '^'^Mgas "^^ Ct.mo't o"Afg„3, 
where Ct,m is the correlation coeffici ent and a the sca tter 
values of temperature and gas mass (jStanek et al.ll2010l ). 

In order to look at the behaviour of such residuals in 
our simulations, we plot in Fig. [3]the residuals in InTmw ver- 
sus the residuals in In A/gas determined at fixed total mass 
Afsoo. The residuals are computed with respect to the best- 
fitting relations computed without imposing the self-similar 
slope. Evaluating the Pearson's correlation coefficient rp for 
the residuals we found them to be uncorrelated. We find 
rp = 0.24 for NR-SV distribution (left panel) and rp = 0.11 
for the CSF-M-W one (right panel). The presence of a weak 
positive correlation explains why the scatter on Yx is, if any, 
slightly larger than that on Afgaa . In order to reproduce m ore 
closely the procedure followed bv lKravtsov et al] (|2006l ). we 
then repeated the analysis by computing the residuals with 
respect the best-fitting self-similar relation. The results in 
this case are shown in Fig. [5] As for the non-radiative runs 
(left panel) we find that A In Tmw and A In A/gas are again 
consistent with being uncorrelated for the NR-SV simula- 
tions, with rp = 0.14. However, such scatters show a more 
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Figure 4. The residuals in \nTmw versus the residuals in InMgas at fixed values of A/500 fof clusters at 2 = 0. Residuals are evaluated 
with respect to the best-fitting M500 — Tmw and M500 — Mgas relations at 2 = 0. Left and right panels are respectively for non-radiative 
(NR-SV) and radiative (CSF-M-W) simulations of clusters of Set 1. Different colours of the symbols correspond to different ranges in 
M500. Cyan triangles: A/500 > 5 X IO^^A/q; Red squares: lO^" < A/500 < 5 X 10^* Mq; Green circles: 5 X 10^^ < A/500 < IO^^'A/q. 



evident sign of anti-correlation for the CSF-M-W simula- 
tions, rp = —0.53. 

In order to make more clear the reason for the anti- 
correlation in the residuals computed with respect to the 
self-similar best-fitting relation for the radiative runs, we 
show with different symbols in Figs. |4] and [5] the points 
corresponding to clusters within different mass ranges: the 
cyan triangles show the most massive objects with M500 > 
5 X 10^'* M0, the green points refer to groups with M500 < 
10^'' Mq, while red squares represent the intermediate clus- 
ters. The only panel showing anti-correlation, the right 
panel of Fig. [S] also shows a clear trend with cluster mass. 
Anti-correlation for CSF-M-W clusters is actually a direct re- 
sult of the offset of scaling relations in temperature and gas 
meiss with respect to the expected self-similar ones. As al- 
ready discussed, in the radiative runs the two relations show 
an opposite behaviour with respect to self-similarity: M- 
Mgas and Af-Tmw are corresponding shallower and steeper 
than the self-similar scaling. This means that massive clus- 
ters in CSF-M-W simulations do have higher gas masses and 
lower temperatures. On the other side, in low mass clusters 
the amount of gas is lower and temperatures are higher than 
expected. Therefore, in this case the presence of an anti- 
correlation is only a spurious effect of imposing the self- 
similar slopes to the best-fitting scaling relations used to 
compute such residuals. Clea r ly, in comparing our results 
with those by iKravtsov et al. I l|2006l ) one should remember 
that their analysis attempted to include observational ef- 
fects in the estimate of the temperature and of the gas mass. 
The analysis presented here, instead, is aimed at quantify- 
ing the intrisic performances of the different mass proxies, 
i.e. before convolving with the effect of measuring temper- 
ature from X-ray spectra and gas mass from an projected 
map of X-ray surface brightness. Quite likely, including such 
observational effects could affect both scaling relations and 



the intrisinc scatter around them through the distribution 
of the residuals. We will address in detail these issues in a 
forthcoming paper (Rasia et al., in preparation). 

In summary, the results presented in this section have 
shown that: (i) M^as has an intrinsic scatter against cluster 
mass slightly smaller than Ix, while both such proxies have 
a significantly lower scatter than Tmw; (ii) Yk is the proxy 
whose scaling against mass has the weakest dependence on 
the inclusion of radiative physics and whose slope is closest 
to the prediction of the self-similar model. The latter result 
is in line with the expectation that Yx is a measure of the 
total thermal energy content of the ICM, which is mostly 
determined by gravitational processes once central cluster 
regions are excluded from the analysis. Moreover, we found 
that residuals in Tmw and Afgas with respect to the best-fit 
relations are uncorrelated independently of the adopted sim- 
ulation, while anti-correlation seen in radiative simulation 
is as a result of the offset of M500 — Tmw and M5Q0 — M^as 
relations with respect to the expected self-similar ones. 

3.3 Stability against changing ICM physics 

In order to further investigate the sensitivity of the differ- 
ent mass proxies on the physics included in the simulation, 
in this section we use results from the analysis of the Set 
2 simulated clusters. As already discussed in Sect. 2, this 
set is characterized by a much smaller number of simulated 
clusters than Set 1, but includes 7 different prescriptions to 
describe the baryonic physics. 

As first test we verify that simulated clusters from Set 1 
and Set 2 do produce consistent results when simulated with 
the same physics (despite the fact they are based on some- 
what different cosmological models) . We show in Fig. [S] re- 
sults from the non-radiative runs (NR-SV and NR-RV) of 
Set 2 and compare them with the best-fitting relation for 



© 0000 RAS, MNRAS 000, 000-000 



Mass proxies from simulations 11 



0.3 



0.2 



0.1 - 



- 



-0.1 



-0.2 - 



-0.3 



1 1 

NR-SV z=9 


1 1 


fixed self-similar slope 

1 1 1 1 1 



-0.3 -0.2 -0.1 0.1 0.2 0.3 

A In Mg3^ 



0.3 



0.2 



0.1 



-0.1 



-0.2 



CSF-M-W z=0 



It • 
-I • 

fixed self-similar slope 



-0.3 

-0.3 -0.2 -0.1 0.1 0.2 0.3 

A In Mgg^ 



Figure 5. The same as in Fig. |4]but for the deviations with respect to the best-fitting relation obtained by fixing the slopes to the 
self-similar values. Colours of the symbols have the same meaning as in Fig. (4] 



the NR-SV version of Set 1 (continuous line). In performing 
this comparison, we should account for the different values of 
the baryon fraction assumed in the two simulation sets (see 
Sect. 2.1). To account for this, we scaled Mgas and Ix for the 
Set 1 simulation NR-SV by multiplying them by the ratio of 
the two baryon fractions: fb,Seti/ fb,Set2 — = 1.2. From this 
plot we note that results for the NR-SV versions of Set 1 and 
Set 2 agree very well with each other. Furthermore, we also 
note that the best fitting relations for the NR-SV and NR-RV 
simulations of clusters of Set 2 are also extremely similar. 
As shown bv lDolag et all (|2005l ), using the time-dependent 
reduced viscosity scheme has the effect of increasin g turbu- 
lent motions in the ICM (see also I Valdarninil [20101 ) . There- 
fore, the stability of results against variations of the adopted 
viscosity scheme (still within SPH hydrodynamics) implies 
that non-thermal pressure support associated to turbulent 
motions have only a very modest impact on the overall gas 
mass and thermal content in our simulated clusters. 

As for radiative simulations, we verified the effect of 
changing feedback efficiency and of including thermal con- 
duction. All such processes are expected to change the gas 
distribution and the temperature structure of the ICM and, 
therefore, to potentially impact on the scaling relations of 
mass proxies against total cluster mass. We plot in Fig. [7] 
such scaling relations obtained the different radiative runs, 
CSF, CSF-C, CSF-M-NW, CSF-M-W and CSF-M-AGN along 
with their best-fitting relations. 

The relations involving gas mass and temperature (left 
and central panel, respectively) show significant dependen- 
cies on the physical processes included in the simulations, 
especially for relatively low mass systems. In general, while 
the effect of including thermal conduction is quite small, the 
separation between simulations with and without AGN feed- 
back is more evident. Indeed, AGN feedback has the effect 
of removing a significant amoun t of gas from the innermost 
regions fe.g.. lFabian et al]|2010l ). thus inducing a decrease of 
Mgas, which is more pronounced for galaxy groups than for 



rich clusters. This effect can be observed in the central panel 
of Fig. [T] where the slope of the best fit relation for AGN 
simulations is of a = 0.81, thus significantly flatter than the 
a — 1 self-similar scaling. The effect of AGN is not only 
to remove a big amount of gas from the center of galaxy 
groups but also to heat up the ICM in the central regions of 
such systems. This can be seen in the M500 — Tmw relation 
shown in the central panel, where the resulting slope of the 
relation is a = 1.73, again quite different from the a = 1.5 
self-similar value. Simulations including the effect of galac- 
tic winds, namely CSF, CSF-C and CSF-M-W, do have a 
slope closer to the self-similar one for both the A/500 — Mgas 
(q ~ 0.9) and the M500 - T^w (a ^ 1.5-1.6) relation. The 
clusters simulated without feedback have an intermediate 
behaviour for both relations. Similarly to what already seen 
for simulations of Set 1, M500 — Tmw is the scaling relation 
characterized by the largest scatter. 

On the other hand, the M500 — Yx relation (right panel 
in Fig. O is confirmed to be the one with the weakest sen- 
sitivity on the baryon physics included in the simulation. 
To further emphasize the different sensitivity that different 
mass proxies have on the ICM physics, we plot in Fig. [S] 
the logarithm of the normalization, log C and the slope a 
of all the best fit scaling relations for the seven different 
physical models simulated for the clusters of Set 2. Tmw and 
Mgas show the largest discrepancies in log C varying from 
14.2 to 14.5 in both relation. The Yx proxy has instead a 

more stable behaviour with values around 14.25 14.35. 

Radiative runs show also a consistent off-set from the self- 
similar slopes in the case of A/gas and Tmw while the slope 
of M500 — Yx relation is in good agreement with the self- 
similar value of a = 0.6 for all the models considered. A 
maximum deviation of ~ 11% is only present in the case of 
CSF-M-AGN feedback. 

In summary, in this section we have shown that the 
Yx mass proxy is the most robust one against changing 
the physical processes which determine the thermodynami- 
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cal properties of the ICM and its scaling with total cluster 
mass has a slope which is always very close to that predicted 
by the self-similar model. 



3.4 Evolution of scaling relations 

For mass proxies to be used in cosmological applications of 
galaxy clusters, one has to precisely calibrate not only nor- 
malization, slope and scatter of their scaling relations with 
total mass for low-redshift systems, but also their evolu- 
tion with redshift. Redshift evolution is especially impor- 
tant when considering distant clusters, located at the high- 
est redshifts, 1, where they have been secured to date 
in statistically complete X-ray surveys. Indeed, as we ap- 
proach the epoch of cluster assembly one expects clusters to 
be characterized by major mergers, which may significantly 
impact on the shaping of such scaling relations. This issue 
is of crucial importance if we want to constrain cosmologi- 
cal parameters by using upcoming and future large cluster 
surveys, based on the next generation of wide field X-ra y 
telescopes (e.g.. lBorgani et al]|2010l : [Cappelluti et al.ir2010t ). 

In this Section we will address the evolution of scal- 
ing relations and of their intrinsic scatter. For an accurate 
study of the evolution of scaling relations we resort to the 
large statistical database of clusters for Set 1 simulations. 
All simulated clusters are analysed at five redshifts: z = G, 
0.25, 0.5, 0.8 and 1. After measuring the normalization C 
and slope a of the scaling relations at each redshift, their 
redshift evolution is described through the relations: 

logC(z) = logCo(^)-h/3i(l + ^); 

a{z) = ao + /32 (1 + z) . (5) 

Similarly to what done by IShort et all l|2010l ). who studied 
the evolution of scaling relations from Millenium simula- 
tions, we take the redshift dependence of Co(z) to be that 
predicted by the self-similar model. Therefore, we expect 
/3i = for a self-similar evolution of the normalization. 

We list in Table U the values of the parameters which 
describe through eqs[5]the evolution of the scaling relations 
between the three mass proxies Afgas , Tmw and Vx with the 
cluster total mass. Moreover, we plot in Fig. [5] the redshift 
dependence of normalization (after accounting for the self- 
similar evolution; upper panels) and slope (lower panels) of 
these scaling relations. The Mgas proxy has a very mild neg- 
ative evolution with redshift (/9i < 0) for both radiative 
and non-radiative simulations. The slope is confirmed to be 
very close to the self-similar value for non-radiative simula- 
tions at all redshift, with a negligible evolution of its value, 
while being shallower for the radiative simulations. As for 
the mass-weighted temperature, a slightly positive evolution 
is detected for the CSF simulations, while the evolution is 
very close to self-similar for NR simulations. 

Quite interestingly, the M500 — scaling relation is the 
one having both the weakest sensitivity to the ICM physics 
at all redshifts and the smallest deviations from the self- 
similar predictions on the evolution of the normalization 
and the value of the slope. These results confirm once again 
the robustness of the Fx mass proxy against variations of 
the physics of the ICM over the whole considered redshift 
range. They further confirm that, since this mass proxy pro- 
vides a measure of the ICM total thermal energy content, it 



can be reliably described on the ground of the predictions 
of the self-similar model. Indeed, this model is based on 
the assumption that gravitational mechanisms determines 
the thermal content of the ICM, an assumption which is 
expected not to be violated by radiative and feedback pro- 
cesses, once we exclude core regions of galaxy clusters. 

A further important aspect that characterizes the be- 
haviour of the three analysed mass proxies at different red- 
shifts is the evolution of their intrinsic scatter. We plot the 
evolution of these intrinsic scatters in Fig. [10] for both CSF 
(red circles) and NR (blue square) simulated clusters. The 
intrinsic scatter for the M500 — Mgas relation, shown in the 
left panel of Fig. 1101 is remarkably low for both simula- 
tions, with a slight positive (negative) evolution for the ra- 
diative (non-radiative) simulations. In general, the intrinsic 
scatter in mass, crinAi, for Afgas varies in the narrow range 
4 — 5 per cent at all considered redshifts. CSF clusters in- 
stead show a mild decrease with a constant 4% scatter below 
z < 0.5. Such small values are expected, since gas mass is 
not sensitive to cluster mergers, which are expected to play 
an increasingly important role at higher redshift. The mid- 
dle panel of Fig. [10] shows instead the evolution of intrinsic 
scatter for the Afsoo — Tmw relation. In this case we note 
that the intrinsic scatter is larger and shows clear sign of a 
positive evolution for both radiative and non-radiative sim- 
ulations, with an increase from crinA/S^ 10 per cent at 2 = 
to ^ 15 per cent at z = 1. The increase of the scatter in 
the Afsoo — Tmw relation agrees with the expectation that 
temperature is more sensitive than gas mass to the pres- 
ence of substructures in the ICM, which are expected to be 
more prominent at high redshift, when cluster mergers are 
more frequent. As for the Afsoo — Yx scaling relation, it has 
a positive evolution which is driven by the positive evolu- 
tion of the scatter in the Afsoo — Tniw relation (right panel 
of Fig. nop . In general the intrinsic scatter increases from 
(Tin A/ — 5 — 6 per cent at low redshift to ;$ 8 per cent at 
z = l. 

The analysis presented in this Section confirms that gas 
mass is the mass proxy behaving best as for the intrinsic 
scatter in its scaling relation with total mass, both for its 
value and for the stability with redshift. 



An important caveat not to over interpret the results 
presented here concern the fact that the small values of in- 
trinsic scatter for all mass proxies have been obtained by 
neglecting observational effects in the measurement of mass 
proxies. Such observational effects are in general expected 
to introduce intrinsic scatter. A typical example is provided 
by the presence of substructures which are not resolved in 
realistic observational conditions. The contribution of the 
X-ray emissivity from the high-density, low-temperature 
substructures impacts both on the estimate of the gas mass 
from the X-ray surface brightness pro file and on the spectro- 
scopic measure ment of temperature (jMazzotta et al.ll2004l : 
lVikhlininl200^ ). In this way, observational determinations of 
gas mass and temperature will have a larger scatter in their 
scaling relation with cluster mass than their mass-weighted 
counterpart do. 
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Figure 6. Scaling relations at z = between the total mass and three mass proxies (Yx and A/gas computed inside -R50O1 Tlnw evaluated 
extracting the central O.lS-Rsoo)- The evolution of the relation is scaled accordingly, with 7 = 2/5,0 and 1, respectively. Best fits of 
non-radiative Set 2 simulations are represented by long dashed green and short dashed blue lines (NR-RV and NR-SV respectively), 
compared with the best fit of the Set 1 NR-SV clusters with red continuous line. For this comparison we scaled the results of Mgas and 
Tmw obtained for clusters of Set 1 simulations to the baryon fraction of Set 2. 




Figure 7. Effect of different physics models on scaling relations between total mass and Yx, Mgas and Tmw at 2 = for the radiative 
simulations of Set 1. The evolution of the relation is scaled accordingly, with 7 = 2/5,0 and 1, respectively. For each simulated physics 
we plot cluster data with purple dots (CSF), cyan squares (CSF-C), blue triangles (CSF-M-NW), green diamonds (CSF-M-W) and red 
triangles (CSF-M-AGN). For clarity, data points are overplotted with the best-fitting relation for each of the five physics, using the same 
colour coding of the points. 



4 CONCLUSIONS 



Galaxy clusters are powerful tools for cosmological studies, 
since the evolution of their mass function constrain the nor- 
malization of the power spectrum, the density parameter 
of dark matter and dark energy as well as the dark energy 
equation of state, through the linear growth rate of density 
perturbations. When observed in the X-ray band, their high 
emissivity allows clusters to be detected so far out to high 
redshifts, z ^ 1.5. However, to fully exploit the potential of 
galaxy clusters as tracers of cosmic evolution, it is necessary 
to understand in detail the relation between mass proxies, 



based on easy-to-measure X-ray observables, and total clus- 
ter mass. 

In this paper we focused on the study of the Mgas Tmw 
and Yx mass proxies and the effect that different physical 
mechanisms have on the mass-observable relations. The aim 
of this analysis was to answer through simulations to three 
questions: (1) Which is the mass proxy that is least sensitive 
to the uncertain knowledge of the physical processes deter- 
mining the thermodynamical structure of the ICM? (2) To 
what extent such mass proxies follow the prediction of the 
self-similar model for the shape of scaling relations and their 
redshift evolution? (3) How large is the intrinsic scatter in 
these scaling relations and how does it evolve with redshift? 
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Figure 8. Dependence of the best-fitting slope a (lower panels) and normalization C (upper panels) of the scaling relations of the three 
mass proxies against total cluster mass. Results for for M500 — Mgas, M500 — Tmw and Afsoo — are shown in the left, central, and 
right panel, respectively. In the lower panel, the horizontal dashed lines mark the values of the slope of the scaling relations predicted 
by the self— similar model. 





log Co 


/3i 


QO 


/32 


CSF-M-W 
NR-SV 


14.307 ± 0.003 
14.152 ±0.001 


A/gas-Msoo 

-0.079 ± 0.015 
-0.086 ± 0.007 


0.927 ±0.003 
0.979 ± 0.002 


0.005 ±0.017 
0.057 ±0.011 


CSF-M-W 
NR-SV 


14.430 ± 0.005 
14.518 ±0.011 


Tmw-Mgoo 

0.059 ±0.025 
0.020 ±0.054 


1.622 ±0.014 
1.524 ±0.012 


-0.047 ± 0.070 
-0.075 ±0.061 


CSF-M-W 
NR-SV 


14.248 ± 0.002 
14.191 ±0.004 


0.007 ± 0.008 
-0.030 ±0.019 


0.591 ±0.002 
0.597 ±0.002 


0.011 ± 0.012 
0.023 ± 0.009 



Table 4. Best fitting values for the parameters determining the evolution of normalization C and slope a of the scaling relations of 
A^gas, Tmw and Yx against Afsoo, as described by eq.l(5j, over the redshift range ^ 2 ^ 1. For each relation the values obtained for 
both CSF-M-W and NR-SV simulations are reported. Reported uncertainties correspond to Icr uncertainties in the linear regression. 



To answer these questions we used galaxy clusters sim- 
ulated with SPH using the GADGET . We combined two 
different sets of simulations: Set 1, composed by more than 
hundred clusters with mass above 5 x lO^/i'^M© was 
used for its high statistics to calibrate scaling relations and 
study their scatter and redshift evolution; Set 2, containing 
many fewer clusters, simulated with seven different different 
physics schemes was instead used to study the effect of: (i) 
thermal conduction, (ii) artificial viscosity, (Hi) cooling and 
star formation, (iv) galactic winds and (v) AGN feedback. 

The main results of our analysis can be summarised as 
follows. 

(1) In non-radiative simulations the relations between 
cluster total mass and the three considered mass proxies 
follow closely the self-similar prediction. 

(2) In radiative simulations the M500 — Mgas and 
Afsoo — Tniw scaling relations show an opposite deviation 



from self-similarity. The net result is a sort of compensating 
effect in Msoo — relation, with a nearly self-similar slope. 

(3) The Afsoo — relation is the most stable against 
changing of the physical processes included in the simula- 
tions, with its slope and evolution being always very close 
to the predictions of the self-similar model, which is based 
on the assumption that gravity only drives the ICM thermo- 
dynamics. Indeed, the Yx proxy is by definition a measure 
of the total thermal content of the ICM, which is dominated 
by the gravitational process of gas accretion. 

(4) Msoo — Afgas is found to be the scaling relation with 
the lowest scatter in mass, 4 — 6 per cent. Moreover, the 
amount of scatter in this case is almost constant with red- 
shift, independently of the simulated physics. 

(5) The scatter in the Msoo — ^x relation is slightly larger 
that for the Msoo — Mgas relation; its intrinsic scatter grows 
with redshift as a consequence of the increase of merging 
events and presence of substructure in the ICM. 
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Figure 9. Rcdshift-dependence of the best fitting parameters wliicli define the scahng relations between mass proxies and Mtot,500- Left 
panel: scaling with Yx,50o; central panel: scaling with A/gas, 500; right panel: scaling with mass-weighted temperature, Tmw (see eq[5ll. 
Results for CSF-M-W and NR-SV simulations are plotted with red circles and blue squares, respectively. For each point we also plot 
the errorbars corresponding to Icr uncertainty from the log-log linear regression, whose size is smaller than the size of the points. The 
horizontal dotted lines in lower panels mark the values of the slope a predicted by the self— similar scaling. As for the amplitude C, its 
redshift dependence mark the degree of deviation from the evolution predicted by the self— similar model. 
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Figure 10. Redshift evolution of the intrinsic scatter ui^m in the scaling relations for simulated clusters of Set 1: results are plotted 
with red circles for CSF-M-W simulations and blue squares for NR-SV simulations. Left, middle and right panels show results for the 
Aftot,500 — Mgas, Mtot,500 — Yx , and Mtot,500 — Tmw, respectively. Note that the y-axis range in the middle panel is twice as large as 
for the other two panels. 



As already discussed in previous sections, the analy- 
sis discus sed in this pap e r diffe rs in spirit from that pre- 
sented bv lKravtsov et all (|2006l ). who instead presented re- 
sults based on the inclusion of projection effects and spec- 
troscopic estimate of the ICM temperature. In order to 
have a first assessment on the effect of carrying out a more 
observationally-oriented analysis, we also computed the Yx 
mass proxy by using the spectroscopic-like definition of tem- 
perature, Tsi. To this purpo se, we followed the procedure de- 
scribed by IVikhlininl l|2006l ). whi ch generalizes the ana lytic 
formula originally introduced by iMazzotta et al] (|2004l ) . to 
include relatively cold clusters with temperature below 3 
keV. We plot in Fig. [TT] the scaling relation between to- 

^ We use the algorithm proposed by IVikhlininI l l2006h that 
resorts to precomputed tables of some parameters for the 
observed spectra as a function of the temperature. Tables 



tal mass and lx,500,si ~ Afgas x Tsi for the radiative and 
non-radiative versions of the clusters of Set 1. As expected 
the scatter in this case is larger than when adopting the 
mass-weighted definition of temperature: it increases from 
fin A/ = 0.06 to 0.10 and from ai^M ~ 0.05 to 0.08 for the 
NR-SV and CSF-M-W simulations, respectively. Moreover 
also a slight deviation from self-similarity is observed for 
non-radiative runs, with a — 0.64, while the slope for radia- 
tive runs, a = 0.60, agrees with the self-similar value. In the 
same figure we also compare our result s with those presented 
in literature, both from simulations (|Kravtsov et al.ll2006l ) 

were generated with the code that is publicly available at 
|http : //hea-www . harvard . edu/'^ alexey/mlxT [ The tables were 
created by fixing Njj = 5 X 10'^'^ for the galactic hydrogen column 
density, solar abundances to the values by [Crevesse & Sauval 
(1998'), the Chandra ACIS-S CCD response function and the pho- 
ton energy range to 0.7 — 10 keV. 
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for radiative simulations (CSF-M-Wwith red circles) and non- 
radiative simulations (NR-SVwith blue squares) of clusters of 
Set 1 at 2 = 0, and results presented in the literature from simula- 
tions and observational data. For this comparison, wc resorte d to 
the spectroscopic— like definition of temperature (Mazzotta c t al.l 
l2004l: IVikhlininll2006l ') for the computation of Yx si in our simu- 
lations. The dot— dashed (green) line corr esponds to the be st— fi t 
relation from the radiative simulation by iKravtsov et al.l lj2006h . 
The continuous (black) line is the best fit relation found by 
ISun et al from the analysis of Chandra data of clusters 

and groups. The lo ng-dashed (grey) line is the best-fit found by 
lOkabe et aO 1 I2OIOI) for the relation between Yx computed from 
XMM-Newton data and weak lensing masses from Subaru obser- 
vations. 



and f rom observational data (I Sun et al ]l2009'; Oka be et al.l 
l20ld ). Our results are in a good agreement although some 
discrepancy is seen at the high mass end of clusters and for 
NR-SV galaxy groups. 

Clearly, a more detailed analysis of how observational 
effects will affect the scaling relations measured from our 
simulated cluster sets requires using a dedicated software , 
like X-ray Map Simula tor (X-MAS; e.g. lGardini et al.ll2004l : 
iRasia et al]|2005l .[ 20081 ) , to extract mock images and spectra 
from simulations and to reduce data using the same proce- 
dure followed for observational data (Rasia et al., in prep.. 
Paper II). 



ACKNOWLEDGEMENTS 

We are greatly indebted to Volker Springel for providing 
us with the non-public version of GADGET . We acknowl- 
edge useful discussions with Gus Evrard, Audrey Kravtsov 
and Chris Miller. This work has been partially supported 
by PRIN-MIUR grant "The Cosmic Cycle of Baryons" and 
PRIN-INAF 2009 grant "Towards an Italian network for 
computational cosmology", by ASI-AAE and ASI-COFIS 
grants and by the INFN-PD51 grant. DF acknowledges the 
support by the European Union and Ministry of Higher Ed- 



ucation, Science and Technology of Slovenia. ER acknowl- 
edges the Michigan Society of Fellow. KD acknowledges the 
financial support by the "HPC-Europa Transnational Ac- 
cess program" and the hospitality of CINECA and of the 
Department of Physics of the University of Trieste. KD also 
acknowledges the support by the DFG Priority Programme 
1177 and additional support by the DFG Cluster of Excel- 
lence "Origin and Structure of the Universe". DF, SB and 
ER acknowledge the hospitality of the Sesto Center for As- 
trophysics (SCfA), where part of the work has been carried 
out. Simulations have been carried out at the CINECA Su- 
percomputing Center (Bologna), with CPU time assigned 
thanks to an INAF-CINECA grant and to an agreement 
between CINECA and the University of Trieste. 



REFERENCES 

Ameglio S., Borgani S., Pierpaoli E., Dolag K., Ettori S., 

Morandi A., 2009, MNRAS, 394, 479 
Arnaud M., Pratt G. W., Piffaretti R., Boehringer H., Cros- 

ton J. H., Pointecouteau E., 2009, ArXiv e-prints 
Balsara D. S., 1995, Journal of Computational Physics, 121, 

357 

Bartelmann M., Steinmetz M., 1996, MNRAS, 283, 431 
Borgani S., Kravtsov A., 2009, ArXiv e-prints, ASL, in 
press 

Borgani S., Murante G., Springel V., Diaferio A., Dolag K., 
Moscardini L., Tormen G., Tornatore L., Tozzi P., 2004, 
MNRAS, 348, 1078 

Borgani S., Rosati P., Sartoris B., Tozzi P., Giacconi R., 
2010, ArXiv e-prints 

Borgani S., Rosati P., Tozzi P., Stanford S. A., Eisenhardt 
P. R., Lidman C, Holden B., Delia Ceca R., Norman C, 
Squires G., 2001, ApJ, 561, 13 

Burns J. O., Hallman E. J., Gantner B., Motl P. M., Nor- 
man M. L., 2008, ApJ, 675, 1125 

Cappelluti N., Predehl P., Boehringer H., Brunner H., 
Brusa M., Burwitz V., et al. 2010, ArXiv e-prints 

Croston J. H., Pratt G. W., Bohringer H., Arnaud M., 
Pointecouteau E., Ponman T. J., Sanderson A. J. R., Tem- 
ple R. F., Bower R. G., Donahue M., 2008, A&A, 487, 431 

Dolag K., Borgani S., Murante G., Springel V., 2009, MN- 
RAS, 399, 497 

Dolag K., Jubelgas M., Springel V., Borgani S., Rasia E., 
2004, ApJL, 606, L97 

Dolag K., Vazza F., Brunetti G., Tormen G., 2005, MN- 
RAS, 364, 753 

Eke V. R., Cole S., Frenk C. S., 1996, MNRAS, 282, 263 

Ettori S., Dolag K., Borgani S., Murante G., 2006, MN- 
RAS, 365, 1021 

Ettori S., Tozzi P., Borgani S., Rosati P., 2004, A&A, 417, 
13 

Evrard A. E., Metzler C. A., Navarro J. F., 1996, ApJ, 469, 
494 

Fabjan D., Borgani S., Tornatore L., Saro A., Murante G., 

Dolag K., 2010, MNRAS, 401, 1670 
Gardini A., Rasia E., Mazzotta P., Tormen G., De Grandi 

S., Moscardini L., 2004, MNRAS, 351, 505 
Gottloeber S., Yepes G., Wagner C, SeviUa R., 2006, 

ArXiv Astrophysics e-prints 



© 0000 RAS, MNRAS 000, 000-000 



Mass proxies from simulations 17 



Grcvesse N., Sauval A. J., 1998, Space Science Reviews, 85, 
161 

Haardt F., Madau P., 1996, ApJ, 461, 20 

Hartley W. G., Gazzola L., Pcarcc F. R., Kay S. T., 

Thomas R A., 2008, MNRAS, 386, 2015 
Henry J. P., Arnaud K. A., 1991, ApJ, 372, 410 
Henry J. P., Evrard A. E., Hoekstra H., Babul A., Mahdavi 

A. , 2009, ApJ, 691, 1307 

Ikebc Y., Rciprich T. H., Bohringer H., Tanaka Y., Ki- 
tayama T., 2002, A&A, 383, 773 
Jeltema T. E., Hallman E. J., Burns J. O., Motl P. M., 

2008, ApJ, 681, 167 

Jubelgas M., Springel V., Dolag K., 2004, MNRAS, 351, 
423 

Kaiser N., 1986, MNRAS, 222, 323 

Kay S. T., da Silva A. C, Aghanim N., Blanchard A., Lid- 
die A. R., Pugct J.-L., Sadat R., Thomas P. A., 2007, 

MNRAS, 377, 317 
Komatsu E., Smith K. M., Dunkley J., Bcmictt C. L., Gold 

B. , Hinshaw G., ct al. 2010, ArXiv e-prints 

Kravtsov A. V., Nagai D., Vikhlinin A. A., 2005, ApJ, 625, 
588 

Kravtsov A. V., Vikhlinin A., Nagai D., 2006, ApJ, 650, 
128 

Lau E. T., Kravtsov A. V., Nagai D., 2009, ApJ, 705, 1129 
Leccardi A., Molendi S., 2008, A&A, 486, 359 
Lima M., Hu W., 2005, Phys. Rev. D, 72, 043006 
Loken C., Norman M. L., Nelson E., Burns J., Bryan G. L., 

Motl P., 2002, ApJ, 579, 571 
Mahdavi A., Hoekstra H., Babul A., Henry J. P., 2008, 

MNRAS, 384, 1567 
Mantz A., Allen S. W., Ebeling H., Rapetti D., Drlica- 

Wagner A., 2010, MNRAS, pp 1030-+ 
Markevitch M., 1998, ApJ, 504, 27 
Maughan B. J., 2007, ApJ, 668, 772 

Mazzotta P., Rasia E., Moscardini L., Tormen G., 2004, 

MNRAS, 354, 10 
Monaghan J. J., 1997, Journal of Computational Physics, 

136, 298 

Morris J. P., Monaghan J. J., 1997, Journal of Computar 

tional Physics, 136, 41 
Muanwong O., Kay S. T., Thomas P. A., 2006, ApJ, 649, 

640 

Nagai D., Kravtsov A. V., Vikhlinin A., 2007, ApJ, 668, 1 
Nagai D., Vikhhnin A., Kravtsov A. V., 2007, ApJ, 655, 98 
O'Hara T. B., Mohr J. J., Bialek J. J., Evrard A. E., 2006, 
ApJ, 639, 64 

Okabc N., Zhang Y., Finogucnov A., Takada M., Smith 

G. P., Umctsu K., Futamasc T., 2010, ArXiv c-prints 
Padovani P., Matteucci F., 1993, ApJ, 416, 26 
PierpaoU E., Borgani S., Scott D., White M., 2003, MN- 
RAS, 342, 163 

Piffaretti R., Arnaud M., Pratt G. W., Pointecouteau E., 

Melin J., 2010, ArXiv e-prints 
Poole G. B., Babul A., McCarthy I. G., Fardal M. A., Bild- 

fell C. J., Quinn T., Mahdavi A., 2007, MNRAS, 380, 437 
Pratt G. W., Bohringer H., Croston J. H., Arnaud M., 

Borgani S., Finoguenov A., Temple R. F., 2007, A&A, 

461, 71 

Pratt G. W., Croston J. H., Arnaud M., Bohringer H., 

2009, A&A, 498, 361 

Press W. H., Teukolsky S. A., Vetterling W. T., Flannery 



B. P., 1992, Numerical recipes in C. The art of scientific 
computing 

Puchwein E., Sijacki D., Springel V., 2008, ApJL, 687, L53 
Rasia E., Mazzotta P., Borgani S., Moscardini L., Dolag K., 

Tormen G., Diaferio A., Murante G., 2005, ApJL, 618, LI 
Rasia E., Mazzotta P., Bourdin H., Borgani S., Tornatore 

L., Ettori S., Dolag K., Moscardini L., 2008, ApJ, 674, 

728 

Rasia E., Tormen G., Moscardini L., 2004, MNRAS, 351, 
237 

Reiprich T., Bohringer H., 2002, ApJ, 567, 716 
Roncarelli M., Ettori S., Dolag K., Moscardini L., Borgani 

S., Murante G., 2006, MNRAS, 373, 1339 
Salpctcr E. E., 1955, ApJ, 121, 161 

Schuecker P., Guzzo L., CoUins C. A., Bohringer H., 2002, 

MNRAS, 335, 807 
Shaw L. D., Holder G. P., Dudley J., 2010, ApJ, 716, 281 
Short C. J., Thomas P. A., Young O. E., Pcarcc F. R., 

Jenkins A., Muanwong O., 2010, ArXiv c-prints 
Sijacki D., Springel V., 2006, MNRAS, 366, 397 
Springel V., 2005, MNRAS, 364, 1105 
Springel V., Di Matteo T., Hernquist L., 2005, MNRAS, 

361, 776 

Springel v., Hernquist L., 2003, MNRAS, 339, 289 
Stanek R., Rasia E., Evrard A. E., Pearce F., Gazzola L., 

2010, ApJ, 715, 1508 
Steinmetz M., 1996, MNRAS, 278, 1005 
Sun M., Voit G. M., Donahue M., Jones C, Forman W., 

Vikhlinin A., 2009, ApJ, 693, 1142 
Sutherland R. S., Dopita M. A., 1993, ApJS, 88, 253 
Tormen G., Bouchet F. R., White S. D. M., 1997, MNRAS, 

286, 865 

Tornatore L., Borgani S., Dolag K., Matteucci F., 2007, 

MNRAS, 382, 1050 
Valdarnini R., 2010, ArXiv c-prints 
Vikhlinin A., 2006, ApJ, 640, 710 

Vikhlinin A., Burenin R. A., Ebeling H., Forman W. R., 
Hornstrup A., Jones C, Kravtsov A. V., Murray S. S., 
Nagai D., Quintana H., Voevodkin A., 2009, ApJ, 692, 
1033 

Vikhlinin A., Kravtsov A., Forman W., Jones C, Marke- 
vitch M., Murray S. S., Van Speybroeck L., 2006, ApJ, 
640, 691 

Vikhlinin A., Kravtsov A. V., Burenin R. A., Ebeling H., 

Forman W. R., Hornstrup A., Jones C, Murray S. S., 
Nagai D., Quintana H., Voevodkin A., 2009, ApJ, 692, 
1060 

Vikhlinin A., Voevodkin A., MuUis C. R., VanSpeybroeck 
L., Quintana H., McNamara B. R., Gioia I., Hornstrup 
A., Henry J. P., Forman W. R., Jones C, 2003, ApJ, 590, 
15 

Voevodkin A., Vikhlinin A., 2004, ApJ, 601, 610 
Voit G. M., 2005a, Advances in Space Research, 36, 701 
Voit G. M., 2005b, Reviews of Modern Physics, 77, 207 
Yang H., Rickor P. M., Sutter P. M., 2009, ApJ, 699, 315 
Yoshida N., Sheth R. K., Diaferio A., 2001, MNRAS, 328, 
669 

Zhang Y., Okabe N., Finoguenov A., Smith G. P., Piffaretti 
R., Valdarnini R., Babul A., Evrard A. E., Mazzotta P., 
Sanderson A. J. R., Marrone D. P., 2010, ApJ, 711, 1033 



© 0000 RAS, MNRAS 000, 000-000 



